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Metallic Glasses are prone to fail mechanically via a shear-banding instability. In a remarkable 
paper Johnson and Samwer demonstrated that this failure enjoys a high degree of universality in 
the sense that a large group of metallic glasses appears to possess a yield-strain that decreases 
with temperature following a — J 12 / 3 law up to logarithmic corrections. In this Letter we offer a 
theoretical derivation of this law. We show that our formula fits very well simulational data on 
typical amorphous solids. 



Introduction: A satisfactory derivation of the rheol- 
ogy of amorphous solids like metallic glasses under vari- 
ous mechanical and magnetic external strains is still far 
from being accomplished. Examples of universal phe- 
nomena and universal relations are rare and far between 
[1]. One of these rare examples is the Johnson-Samwer 
T 2 / 3 law [2] which pertains to the temperature depen- 
dence of the yield-strain (the value of the strain where 
the material fails via plastic instabilities). Examining a 
large group of metallic glasses these authors proposed the 
law 

7y (T, 7 ) = 7y (T = 0,7 = 0) {l - [ATln( Wo /C 7 )] 2/3 } , 

(1) 

where 7 Y (T = 0,7 = 0) is the yield-stress at athermal 
quasi-static conditions, 7 is the external strain rate, A 
and C are material constants and wo is a microscopic in- 
verse time scale. Johnson and Samwer offered a deriva- 
tion of this law [2] based on Frenkel's theory [3] which is 
appropriate for an infinite crystals having a periodic elas- 
tic energy density, leaving a derivation which is proper 
for an amorphous solid for the future. This task was 
picked up in Rcf. [4]. These authors recognized that in 
amorphous solids plastic yielding follows a saddle-node 
bifurcation where an eigenvalue Xp of the system's Hes- 
sian hits zero at a value of the strain 7 = 7 p. Such 
a bifurcation leads to a local rearrangement which is 
characterized by an energy barrier SE which scales like 
SE ~ (7P — 7) 3 / 2 . Assuming that this energy barrier 
is much larger than kT (k being Boltzmann's constant), 
Ref. [4] showed that a temperature induced barrier cross- 
ing would result in a law of the form of Eq. 1. In this 
Letter we argue that the energy barrier that needs to be 
crossed to achieve shear localization is sub-extensive; on 
the other hand the energy barrier associated with a single 
localized plastic instability is minute, and in the thermo- 
dynamic limit it scales like 1/N where N is the number 
of particles [5]. Therefore the conditions at which the 
system yields via a coherent shear localized band need to 
be reconsidered. 

Yielding via shear localization: In recent work it 
was argued that when the external strain approaches the 
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FIG. 1: (Color Online). The total plastic energy Eq. (3) 
for the creation of an array of quadrupoles with density p for 
three values of 7: 7 = 7y — 0.1 (upper curve), 7 = 7 Y — 0.05 
(middle curve), and 7 = 7 Y (lower curve). In the present 
case 7 Y = 0.07. To generate this picture we use the measured 
constants £ « 37.2, v m 0.31, e* m 0.082 and a = 1.83. 
Finally U p m 0.22. 

shear-localization yield-strain the nature of plastic insta- 
bilities can change qualitatively [6, 7]. At low external 
strains a plastic instability results in a local rearrange- 
ment such that the non-affinc displacement of particles 
can be very well modeled by the displacement field asso- 
ciated with a single Eshelby inclusion in an elastic matrix 
[8] . In 2D this field assumes a quadrupolar structure. At 
larger values of the strain the instability results in 2D in 
a highly correlated array of Eshelby quadrupoles that are 
aligned at 45 degrees to the principal stress axis. All the 
quadrupoles are in phase and in total they result in shear 
localization in a narrow strip. One could show that this 
highly correlated array is a minimum energy state which 
depends on the density p of the quadrupoles. For 7 < 7 Y 
the only solution is p = 0, i.e. isolated qudrupolcs, but 
at 7 = 7 Y a bifurcation opens up a new solution with a 
finite density, cf. Fig. 1. Under conditions of athermal 
quasistatic straining (AQS), it was shown that the yield 
strain 7 Y is given by the expression [7] 

7y = 7 Y (^ = 0,7 = 0) ^2(T~^)- (2) 
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where e* is the cigenstrain induced by single plastic insta- 
bility in the background elastic matrix whose Poisson's 
ratio is v. Each quadrupole can be very well modeled as 
an Eshelby inclusion with eigenstrain e* and core size a 
[8] . The energy density cost of creating a linear array of 
Af quadrupoles all with the same orientation, separated 
by distance R = L/J\f (in a 2-dimensional system of size 
L 2 ) was computed analytically [7] in the form 
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where U p = [£ 7r(e*) 2 ]/[4(l - v 2 )\ with £ being Young's 
modulus; while B = 4((2) and C = 6C(4), where the 
Riemann zeta functions are respectively £(2) = 7r 2 /6 and 
C(4) = 7r 4 /90. This energy is shown in Fig. 1. Examining 
this energy we realize that in AQS conditions yield via 
shear localization can occur only at 7 = 7 Y . With finite 
temperature we can have thermally-assisted transitions 
which we consider next. 

Thermally assisted plastic yield: Considering the 
yield stress at any finite temperature T / 0; we real- 
ize that thermal fluctuations can always (if given enough 
time) cross any energy barrier and therefore under con- 
ditions of quasistatic straining the yield stress should al- 
ways vanish 



7y (T, 7 = 0) = 0. 



(4) 



Accordingly, for strictly quasistatic straining 7 = the 
limit T — > yield-stress is different from the T = yield- 
stress. This jump in the limit was observed in [2] and 
was interpreted as a "quantum effect" . There is no need 
to invoke quantum mechanics to understand this simple 
issue. 

In reality one always strains at some finite value of 7. 
Then thermal fluctuations can lead to barrier crossing if 
the timescale for such crossing r is smaller or similar to 
the straining time scale 7 , leading to the condition for 
yielding due to thermal fluctuations given by 



7-7 = O(l). 



(5) 



If we had a process in which every particle in the system 
acted independently, it would be enough to compute the 
energy barrier per particles, say Ub, and compare it to 
the thermal energy, to be used in the Arrhenius form 

r = To exp [Ub/kT], for independent particle motion 

(6) 

where tq is the inverse of the attempt frequency for bar- 
rier crossing. In fact, the situation here is more complex. 
The way that shear localization occurs in practice is that 
there exists a mode that gets localized on N* <C N parti- 
cles which are involved in the first quadrupole that forms 
and then the rest of the linear array of qudrupoles forms 
instantly, as seen in Fig. 2. Thus the rate determin- 
ing step is the creation of the first quadrupole, and we 
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FIG. 2: (Color Online). Visualization of the process of shear 
localization in the AQS simulations whose full details can be 
found in Ref. [7]. Even with the quasistatic protocol with 
arbitrary long waiting times if necessary one cannot resolve 
the formation of the full structure of shear localized band 
of quadrupoles (right panel) from the creation of a single 
quadrupole (left panel). The right panel appears instantly 
after the left panel. 



estimate the appropriate times scale as 

r = To exp [N*Ub/kT], for concerted particle motion 

(7) 

Combing these results gives the expression for yielding 

as 



To7exp[Ar*CV£;T] = O(l). 



(8) 



Note that TV* remains independent of N in the ther- 
modynamic limit, representing the number of particles 
involved in the concerted barrier crossing of the creation 
of one quadrupolar non-affine displacement field. We ex- 
pect N* to be of the order of 100, give or take a factor 
of 2. 

We now use these result to estimate how finite tem- 
peratures and finite strain rates change the value of the 
AQS yield strain due to shear localization. 

Calculating the energy barrier: Examining Eq. 
(3) we see that the plastic energy density due to shear 
localization at strain 7 can be written in terms of the di- 
mensionless variable x = pa. The height of the barrier per 
particle Ub is determined by the energy density E* / La at 
x* = p*a which can be found from d(E/ La)/dx = with 
the result that the barrier occurs at 



x* =p*a= v/[l- 7 / 7y ]/2/7r, (9) 
and its strip energy density is given by 

V2U P 



E* /La 
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(10) 



Eq. (10) implies that the barrier energy per particle in 
the strip is 

U b = {E* I La) In = [V2iy (3™)][1 - 7 /7.] (3/2) , (11) 

where n is the number density of particles in the strip of 
dimension La. 
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Finally we need to estimate To, the inverse attempt 
frequency. There are two candidates. The first is the 
eigenvalue of the Hessian matrix A p which is associated 
with the eigenfunction that gets localized. The square 
root of A p is an inverse time scale. This time scale how- 
ever diverges to infinity when 7 is close to 7 Y , and the 
purely thermal time scale becomes shorter and more rel- 
evant. We thus argue here that the relevant time scale 
can be estimated from the thermal fluctuations of the in- 
dividual particles in the solid which determines the time 
scale of a concerted motion of N* particles. For a single 
particle the thermal time scale is l/v where I « 
is the typical distance between particles and the typical 
velocity v can be found from the equipartition theorem 
m(v 2 } = kT. Thus for N* particles moving together we 
estimate 



to 



N*yJm/(nkT). 



(12) 



Now combining Eq. (8) for the yield strain with Eq. (11) 
for the barrier height and Eq. 12 for the bare time scale 
we finally find 

2/3 

Xttti fa' I ' 



3imkT 



V2N*U P 



log 



2/3 



yj nkT Jr 



(13) 

We note that at this point only N* is not known with 
certainty, only in order of magnitude. The expression 
has the T^ 2 / 3 ) temperature reduction in the yield strain 
found experimentally by Johnson and Samwer, cf. Eq. 
(1), including the logarithmic correction term. 

Comparison with numerical simulations: We 
have performed 2D Molecular Dynamics simulations on 
a binary system which is an excellent glass former and is 
known to have a quasi-crystalline ground state. Each 
atom in the system is labeled as either "small" (S) or 
"large" (L) and all the particles interact via Lennard 
Jones (LJ) potential. All distances are normalized by 
(Jsl, the distance at which the LJ potential between 
the two species becomes zero and the energy is normal- 
ized by esL which is the interaction energy between two 
species. For detailed information on the model poten- 
tial and its properties, we refer the reader to Ref [10]. 
The number of particles taken in all our simulations is 
100489 at a number density n = 0.985 with a particle 
ratio Nl/N s = (1 + y/E}/4. The mode coupling temper- 
ature Tmct for this system is known to reside close to 
0.325esL/fcs. All particles have identical mass and 
hence the time is normalized to to = \J csl&sl 2 /mo- 
For the sake of computational efficiency, the interaction 
potential is smoothly truncated to zero along with its 
first two derivatives at a cut-off distance r c = 2.5<jsl- 
To prepare the glasses, we start from a well equilibrated 
liquids at a high temperature of 1.2esi,/kB which are su- 
percooled to 0.35esi/fcs at a reduced quenching rate of 
3.4 x 10 -3 ig . We then equilibrate these supercooled liq- 
uids for times greater than 10r re i, where r rc i is the time 
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FIG. 3: (Color Online). Typical stress vs. strain curve ob- 
tained at finite temperature and finite strain rate. The yield 
stress was estimated by fitting a cubic to the curve (see inset) 
and finding the maximum of the curve. 



taken for the self intermediate scattering function to be- 
come 1% of its initial value. Following this equilibration, 
we quench these supercooled liquids deep into the glassy 
regime at a temperature of O.OSesL/ks at a much slower 
quenching rate of 3.2 x 10 -6 ijj" . 

We perform simple shear loading experiments by inte- 
grating the SLLOD system of equations [11] along with 
Lees-Edwards periodic boundary conditions [12]. As can 
be expected, work done on the system will lead to heat 
dissipation and hence we keep our system connected to 
the thermostat in order to keep the temperature constant 
throughout the mechanical deformation. We employ a 
strain rate 7 = 10~ 5 in all our loading experiments. A 
typical stress vs. strain curve is shown in Fig. 3. The 
data of stress vs strain were fitted to a cubic in 7 and 
the yield strain was estimated from the maximum of the 
curve. 

To compare with the theory above we need also "fy at 
zero temperature and the value of U p . We obtained these 
by minimizing the energy of one of the glasses generated 
as explained above to T = 0. Performing athcrmal quasi- 
static simulations on this sample wc determined all the 
parameters appearing in Eq.(13), with the result "fy ~ 
0.06 and U p = 0.22. Armed with these and all the other 
known numbers we compare the data for the temperature 
dependent yield strains to the prediction of Eq. (13) 
using N* to get a best fit, which is obtained for N* ss 250. 
The fit is shown in Fig. 4. Having in mind the inherent 
unknown factors hidden in the estimates Eqs. (5), (6) and 
(12), we find the fit very satisfactory. 

In summary, we have used the analytically computed 
energy associated with a correlated series of quadrupolar 
structures that add up to a shear localizing instability 
to study the thermally assisted yield strain associated 
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FIG. 4: (Color Online). Comparison of the simulational data 
of 7(T, 7) to the prediction Eq. (13). The dotted line was 
obtained by finding the best fit for N* . The error bars indicate 
averaging over three independent realizations. 

with plastic failure in two dimensions. This energy is 
sub-extensive and its strip-density remains unchanged in 
the thermodynamic limit. We could derive the Johnson- 
Samwer T 2//3 law essentially without a free parameter ex- 
cept for an uncertainty regarding the value of ./V* . The 
fitted value of this number appears to be in the right or- 
der of magnitude, lending strong support to the approach 
detailed above. 
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